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Using the Markovian method, we study the stochastic nature of electrical discharge current fluc- 
tuations in plasma. Sinusoidal trends are extracted from the data set by the Fourier-Detrended 
Fluctuation analysis and consequently cleaned data is retrieved. We determine the Markov time 
scale of the detrended data set and show that it is almost a monotonic increasing function of dis- 
charge current intensity. We also estimate the Kramers-Moyal's coefficients of the discharge current 
and derive the corresponding Fokker-Planck equation. The obtained Langevin equation enables us 
to reconstruct discharge time series with similar statistical properties compared with the observed 
in the experiment. Finally the multifractal behavior of reconstructed time series using its Keramers- 
Moyal's coefficients and original data set are investigated. 



I. INTRODUCTION 



Many natural phenomena are identified by a degree 
of stochasticity. Turbulent flows, seismic recordings and 
plasma fluid are but a few examples of such phenomena 
[1-28]. Recently a robust statistical method has been 
developed to explore an effective equation that can re- 
produce stochastic data with an accuracy comparable to 
the measured one [18,19,29-34]. As in many early re- 
searches has been confirmed, one may utilize it to: 

1 : reconstruct the original process with similar statisti- 
cal properties, and 

2 : understand the nature and properties of the stochastic 
process [18-22,29,30]. 

Interpretation and estimation of physical and chemi- 
cal properties of plasma fluid have been one of the main 
research areas in the science of electromagnetic hydro- 
dynamics. Fluctuations of electric and magnetic fields 
of plasma, spectral density, logistic mapping and nonlin- 
earity of ionization wave have been investigated in Refs. 
[10-17,35-37]. It is well-known that discharge fluctua- 
tions in the plasma often exhibits irregular and complex 
behavior. 

As we discussed in Ref. [17], because of some limi- 
tations in the experimental setup and data acquisition 
methods the measured plasma fluctuations may be af- 
fected by some trends such as; alternative power current 
oscillation, noise due to the electronic instruments and 
also fluctuations of striation areas near the anode and 
cathode plates. Consequently, inferring the reliable sta- 
tistical properties of original fluctuations and to avoid 
spurious detection of correlations need to apply robust 
methods. Data set used through this paper were collected 
by a positive column of a hot cathode discharge in he- 
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FIG. 1. Typical detrended electrical discharge current fluc- 
tuations in the plasma as a function of time. 



lium plasma. Details of our experimental setup is given in 
Ref. [17]. Cleaned data were constructed by applying the 
Fourier-Detrended Fluctuation analysis. Figure 1 shows 
typical detrended discharge current fluctuation. The size 
of recorded data points in each sampling is about 10^. 

Here using the Makovian method, we explore the sta- 
tistical properties of discharge current fluctuations. We 
show that this approach will profoundly gives deep in- 
sight through the electromagnetic and hydrodynamics of 
discharge current in the ionized fluid. By addressing the 
implications dictated in [18,19] a Fokker-Planck evolution 
operator and Langevin equation will be found. 

The rest of this paper is organized as follows: Section 
II is devoted to a brief summary of the most important 
notions and theorems on Markovian method and their ap- 
plication to the analysis of empirical data. Using the like- 
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FIG. 2. Normalized correlation function for plasma cleaned 
data for discharge current with / — 50mA (left panel) and for 
completely random data (right panel) 

lihood statistics, wc determined Markov time scale. Sec- 
tion III contains the main results of our analysis and es- 
timate the Fokker-Planck and Langevin equations which 
govern the probability density function and stochastic 
variable (discharge current), respectively. Comparison 
between statistical and multifractal properties of original 
and reconstructed time series are also given in Section III. 
Section IV closes with a discussion and conclusion of the 
present results. 



II. MARKOVIAN NATURE OF DATA SET 

Let us begin with one of the simplest and most impor- 
tant statistical quantity of a given time series, that is the 
correlation function, which is written as 



C(Xi(ti),X2(i2)) = {Xiih)X2{t2)) 



(1) 



here the sign (.), shows the ensemble averaging. We have 
set the average of time series equal to zero. For stationary 
time series the correlation function depends on only the 
separation time scale, which means that 



C(r) =C(Xi(ti),X2(t2)) 

= {X,{h)X2(tl+T)). 



(2) 



In the presence of any trends and nonstationarity, cor- 
relation function depends not only to the time separation 
(t2 — ti), but also to the starting and finishing times, 
namely ti and t2, respectively. As demonstrated in Ref 
[17], the underlying detrended data for discharge current 
behaves as a stationary signal. The correlation function 
for plasma detrended data with / = 50mA is plotted in 
the left hand side of the Figure 2. This figure confirms the 
underlying data sets behave as an anti-correlated series 
which has been observed in [17] using another method. 
We also present the same plot for a pure random signal, 
i.e. with Hurst exponent H ~ 0.5, in the right hand side 
of the Figure 2 for comparison. 

As mentioned in introduction, here we use the Marko- 
vian method to explore the nature of stochasticity of dis- 
charge current fluctuations in plasma. To investigate the 
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FIG. 3. Markov time scale as a function of discharge cur- 
rent intensity. The unit of vertical axis is sec/44100. 



markovian nature of data, we briefly summarize the no- 
tions and theorems which will be importance for our sta- 
tistical analysis of cleaned data set. For further details 
on Markov processes we refer the reader to the references 
[33,38-42]. 

Represent the discharge current fluctuations as a func- 
tion of time by X{t) and define x{t) = X{t)/a, where 
a is the standard deviation of discharge current fluctu- 
ations. Fundamental quantities related to the Markov 
processes are conditional probability density functions. 
The conditional probability density function (CPDF), 
p(.T2, t2la;i, ii), is defined as 



p{x2,t2\xi,ti) 



p(x2,t2;xi,ti) 
p{xi,ti) 



(3) 



where p{x2 ,t2',Xi,ti) is the joint probability density func- 
tion ( JPDF) , describing the probability of finding simul- 
taneously, Xi at scale(time), ti, and X2 at scale(time), 
t2. Higher order conditional probability densities can be 
defined in an analogous way 



p{xN,tN\xN-l,tN-i; ...■,Xi,ti) 



p{xN-l,tN-i; Xi,ti) 

(4) 



where p(xjv, ^tv; a^w-i, ijv-i; a;i, ti) is A^-point joint 
probability density function. Intuitively, the physical in- 
terpretation of a Markov process is that it "forgets its 
past," or, in other words, only the most nearby condition- 
ing, namely xjv-i at tjy-i, is relevant to the probability 
of finding a fluctuation xn at tjv. Hence, in the Markov 
process the ability to predict the value of a; at will not be 
enhanced by knowing its values in the steps prior to the 
most recent one. So an important simplification that is 
made for a Markov process is that, the conditional mul- 
tivariate joint PDF is written in terms of the products of 
simple two parameter conditional PDF's [38] as 
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FIG. 4. Upper panel shows contour plots of the conditional 
PDF, p(x3, ts\xi,ti). The solid and dashed line correspond to 
the left and right hand side of Eq. (11) for t-^ —t\ = 2 x tMarkov , 
respectively. Inner contours are a cutting of PDF at 0.08 
level and outer contours correspond to 0.005 level. Lower 
panel corresponds to the cuts through the conditional PDF 
for xi = ±1.25(7. 

p{xN,tN;XN-l,tN-i; • • • ; X2,t2\xi,ti) 
N 

= Y\_PiXi,ti\xi-i,ti^l) (5) 

i=2 

To investigate whether underlying signal is a Markov 
process, one should tests the Eq. (5). But in practice for 
large values of N, is beyond the current computational 
capability. For = 3 (three points or events), however, 
the condition will be 

p{x3,t3\x2,t2;xi,ti) =p(x3,t3|a;2,t2) (6) 

which should hold for any value of t2 in the interval 
ti < t2 < ts- A process is then Markovian if the Eq. (6) 
is satisfied for a certain time separation ^3 — ^2, in which 
case, we define the Markov time scale as iiviarkov =^3 — ^2- 
For simplicity, we let t2 —ti = ^3 — ^2- Thus, to compute 
the ^Markov wc usc a fundamental theory of probability ac- 
cording to which we write any three-point PDF in terms 
of the conditional probability functions as 

p{x3,t3;X2,t2;Xi,ti) 

= Pix3,t3\x2,t2;Xi,ti)p{x2,t2; Xi,ti) (7) 



FIG. 5. Drift, D'^^\x), diffusion and D^'^\x) coefficients 
for data set with / — 50mA. 

Using the properties of Markov processes to substitute 
Eq. (7), we obtain 

PMa.i{x3,h;X2,t2;Xi,ti) 

= Pix3,t3\x2, t2)p{x2,h\ Xi,ti) (8) 

In order to check the condition for the data being a 
Markov process, we must compute the three-point JPDF 
through Eq. (7) and compare the result with Eq. (8). 
The first step in this direction is determining the quality 
of the fit through the least-squared fitting quantity 
defined by [4] 

= J dxidx2dx3[p{x3,t3;x2,t2;xi,ti) 

~PMi,v{x3,t3;X2,t2;Xl,ti)f/ [disjoint + f^Mar] (9) 

where C3_jQi„^ and cr^^^j. are the variances 
ofp{x3,t3;x2,t2;xi,ti) a.ndpMa^r{x3,t3]X2,t2;xi,ti), re- 
spectively. To compute the Markov time scale, we use the 
Likelihood statistical analysis [43]. In the absence of a 
prior constraint, the probability of the set of three-points 
JPDF is given by a product of Gaussian functions 

xi,X2,X3 y ^"'(.''■s-joint + ^MarJ 
[ \p{X3,t3; X2,t2; Xi,ti) - PMar(a:3, X2, ^2! Xi,ti)]'^ \ 

''''^V 2(a|_j„i„,+<^J J 

(10) 

This probability distribution must be normalized. Evi- 
dently, when, for a set of values of the parameters, the 
is minimum, the probability is maximum. The mini- 
mum value of xi (xt — X^/-^i with TV being the number 
of degree of freedom) corresponds to ^Markov for different 
value of electrical discharge current intensities. The val- 
ues of Markov time scales, iMarkov in terms of discharge 
current intensity have been plotted in Figure 3. It must 
be pointed out that, the unit of iMarkov can be changed 
to the units of Seconds, using the rate of digitalization 



TABLE I. The values of Kramers-Moyal coefficients for 
data set at different discharge current intensities. 







D^''\x) 


50mA 


-0.160 X 


0.090 + 0.003 a: + 0.070 x^ 


60mA 


-0.058 X 


0.026 + 0.002 a: + 0.030 x^ 


fOOmA 


-0.052 X 


0.026 + 0.002 a- + 0.026 x^ 


f20mA 


-0.028 X 


0.0f3 + 0.00f a: + 0.0f4 x^ 


f40mA 


-0.0f7a; 


0.008 + 0.00f a + 0.009 x^ 


fSOmA 


-0.0f7a; 


0.008 + 0.00f a + 0.009 x^ 


2f0mA 


-0.0f6 X 


0.009 + 0.00f a + 0.008 x^ 



in the experimental setup, 44100 sample/ sec (for simpli- 
fication we use "sec" as an abbreviation of "second" for 
unit of time in whole of paper). 

One can write Eq. (8) as an integral equation, which is 
well-known as the Chapman-Kolmogorov (CK) equation 

p{X'i,h\xi,ti) = j dx2 P{x3,t3\x2,t2) p{x2,t2\xi,Xi) 

(11) 

We have checked the validity of the CK equation for 
describing the time scale separation of ti and ^2 being 
equal to the Markov time scale. This is shown in Fig- 
ure 4 (for the data set with electrical current intensity, 
/ = 50 mA). In this figure, the upper panel shows the 
contour plot of identification of the left (solid line) and 
right (dashed line) sides of Eq. (11) for two levels, 0.080 
(inner contour) and 0.005 (outer contour). The condi- 
tional PDF p{x3,t3\xi,ti), for xi = ±1.25tT, are shown 
in the lower panel. All the scales are measured in unit 
of the standard deviation of the discharge current fluc- 
tuations. We must point out that if all situations to be 
same as our experimental setup such as pressure, current 
intensity and so on, one can expect that all values de- 
rived by Markov analysis would be repeated. The value 
of Markov time scale increases as discharge current in- 
tensity increases (see Figure 3). It seems that by increas- 
ing the current intensity, charges become more energetic, 
therefore their effective cross-section will decrease and 
hence increasing their memory. 

Up to now we determined the Markov time scale for 
each cleaned data set over which time series behaves as a 
Markov process. In the next section we will turn to the 
deriving master and stochastic equations governing the 
evolution of probability density function and fluctuation 
itself, respectively. 
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FIG. 6. Left panel corresponds to the conditional proba- 
bility density function determined by analytical formula, Eq. 
(17) (solid line) and directly computed by original cleaned 
data (symbol) for / — 50mA. Right panel shows the com- 
parison between the conditional probability density function 
determined by generated data using Eq. (16) (triangle sym- 
bol) and our initial cleaned data (circle symbol). In each 
panel, the plots from left to right correspond to the cut for 
xi = —0.5a, x\ = 0.0 and x\ = +0.5(t level, respectively. To 
make more obvious, we shifted the value of 2:2 for each plot. 
We took T = tMarkov, whcrc tiviarkov IS the Markov time scale 
of data set. 

III. LANGEVIN EQUATION: EVOLUTION 
EQUATION TO DESCRIBE THE PLASMA 
DISCHARGE CURRENT FLUCTUATIONS 

The Markovian nature of the plasma electrical dis- 
charge fluctuations enables us to derive a Fokker-Planck 
equation - a truncated Kramcrs-Moyal equation - for the 
evolution of the PDF p{x,t), in terms of time t. The 
Chapman-Kolmogorov (CK) equation, formulated in dif- 
ferential form, yields the following Kramers-Moyal (KM) 
expansion [38] 

-p(2;,i) = ^[- — j [D(-\x,t)pix,t)] (12) 

where {x,t) are called as the Kramers-Moyal's coeffi- 
cients. These coefficients can be estimated directly from 
the moments, and the conditional probability dis- 

tributions as 

D^'^\x,t) = ^ lim Af(") (13) 

M^") = ^ y dx'{x' - x)>(x', t + Ai|x, t) (14) 

For a general stochastic process, all Kramers-Moyal's 
coefficients are different from zero. According to 
Pawula's theorem, however, the Kramers-Moyal expan- 
sion stops after the second term, provided that the fourth 
order coefficient D^^\x,t) vanishes. In that case, the 
Kramers-Moyal expansion reduces to a Fokker-Planck 
equation (also known as the backwards or second Kol- 
mogorov equation) [38] 



dt 



(15) 



Also the evolution equation for conditional probability 
density function is given by the above equation except 
that p{x,t) is replaced by p{x,t\xi,ti). Here D^^^ is 
known as the drift term and D*-^^ as diffusion term which 
represents the stochastic part. The Fokker-Planck equa- 
tion describes the evolution of probability density func- 
tion of a stochastic process generated by the Langevin 
equation (we use the Ito's definition) [38] 



^x(<) = {x, t) + ^D(^)ix,t)f{t) (16) 

where f{t) is a random force, i.e. (5-correlatcd white 
noise in t with zero mean and gaussian distribution, 
(fit) fit')) = 2(J(t - t'). Using Eqs. (13) and (14), for 
collected data sets, wc calculate drift, D^^\ and diffu- 
sion, D^^\ coefficients, shown in Figure 5. It turns out 
that the drift coefhcient D^^^ is a linear function in t, 
whereas the diffusion coefficient Z)'-^-' is a quadratic func- 
tion. For large values of x, our estimations become poor, 
the uncertainty increases, so we truncate our estimations 
up to 3(7 of fluctuations as indicated in Figure 5. 

The functional feature of drift and diffusion coefficients 
for different electrical discharge data sets are reported in 
Table I. To ensure that Kramers-Moyal expansion (Eq. 
(12)) reduces to a Fokker-Planck equation (Eq. (15)), we 
compute fourth-order coefficient D^^\ In our analysis, 
£)(4) ^ 10~^D(^). One must point out that, however the 
fourth-order Kramers-Moyal's coefhcient is not so small, 
but in the current analysis, this doesn't make measur- 
able uncertainty in our results (see below). Furthermore, 
using Eq. (16), it becomes clear that we are able to sep- 
arate the deterministic and the noisy components of the 
fluctuations in terms of the coefficients Z?'^' and D'--^\ 
According to the values of the Kramers-Moyal's coeffi- 
cients reported in Table I, it is possible to reconstruct 
discharge current fluctuations at arbitrary current inten- 
sity using Eqs. (15) and (16) [18]. 

Now let us have a comparison of the statistical proper- 
ties of reconstructed data using Eq. (16) with the original 
fluctuations. For this purpose, we rely on the solution of 
Fokker-Planck equation for conditional probability func- 
tion (same as Eq. (15) for infinitesimally small step t) 
which is given by [38] 



p{x2,t + T\xi,t)^ 



1 



2^7r£>(2)(x2,t)T 

(x2-Xi-i?W(x2,t)T)2 



X exp 



4L>(2)(x2,0l 



(17) 



Left panel of Figure 6 shows conditional probability 
density function computed by the above equation and 




FIG. 7. Upper panel indicates the structure function ver- 
sus T. Lower left panel shows the scaling exponent of S'q(r) as 
a function of moment for original cleaned plasma fluctuations 
(filled symbol) and solid line corresponds to the fitting formula 
derived by Kramers-Moyal's coefficients for multi-fractal an- 
ti-correlated signal with H — 0.42 (see Eq. (27)). Also in this 
panel, dashed line corresponds to a mono-fractal anti-corre- 
lated series. Lower right panel indicates ({q) versus q. Here 
we chose the data set with / — 50mA. 



directly calculated from the original detrended data set 
for / = 50mA. The plot from left to right correspond to 
xi = —0.5(7, xi = 0.0 and xi ~ +0.5(7 level, respectively. 
We also compute the conditional probability using recon- 
structed fluctuations via Eq. (16) and compare it with 
the same one for original cleaned data at three mentioned 
levels for xi. We took r — tMarkov for all plots in Figure 
6. 

To check the multifractal nature of reconstructed time 
scries, we investigate the Markovian nature of the incre- 
ments which is defined as: Ax{t) = x{t + r) — x{t). Ac- 
cording to the mentioned procedure, we can determine 
the Markov time scales for the increments and calcu- 
late the Kramers-Moyal's coefficients. Likelihood analy- 
sis confirms that, the increment signal of plasma fluctua- 
tions for all electrical current intensities are also Markov 
processes. The Fokker-Planck equation for probability 
density function of the increment is given by [44,45] 

d 



-t^p{Ax,t) = 

OT 



OA 



dAa 



D^'^\Ax,t)\p(Ax,t) (18) 



the negative sign of the left-hand side of Eq. (18) is due 
to the direction of the cascade from large to smaller time 
scales r. The corresponding Langevin equation can be 
read as 



-t^Ax{t) = D<-^^ (Ax, r) + ^i:)(2)(Ax, t)/(t) (19) 

where /(r) is the same as random function in equation 
16. For time series with scahng correlations the Drift 
and diffusion coefficients of increment are formulated as 
[44-46] 

D'^^\Ax,t) ~ -HAx 

£)(2^(Ax,t) ~ 6Ax2 (20) 

Using Eqs. (18) and (20) we obtain the evolution of 
structure functions ( Sq{T) = (|Ax(r)|'') = {\x{t + t) — 
x(t)\'')) as follows 

-t-^{\Ax{t)\i) = q(\Ax{T)\'i-^D<^^\Ax.T)) 

OT 

+q{q - 1){\Ax{t)\'^-''D^^\Ax, t)) (21) 
by substituting the Eqs. (20) in Eq. (21) we find 

r-^{\Ax{r)n = [qH-bqiq~l)]{\AxiT)\'^) (22) 

the above equation implies scaling behavior for moments 
of increments, structure function as 

Sg{r) EE (|Ax(t)|«) = {\x{t + r) - x(t)|«) ~ t^^"^ (23) 

According to Eqs. (22) and (23), the corresponding 
scaling exponent in general case can be read as 

= qH- bq{q - 1) (24) 

For mono- and multi-fractal processes the exponent ^(g) 
have linear and non-linear behavior with g, respectively. 
It must point out that H is nothing except the underly- 
ing fluctuations's Hurst exponent [47-50]. To check the 
estimated scaling exponent ^(g) ( Eq. (24) ) with origi- 
nal time series we use the extended self similarity (ESS) 
method [51,52]. In the ESS method, the log-log plot of 
Sq (t) as a function of specific order of structure function, 
namely S^{t)^ usually shows an extended scaling regime 

Sq{T) ~ 53(r)?(«) (25) 

For any Gaussian process, the exponent in the above 
equation is given by C,{q) = g/3 [51,52]. Any deviation 
from this relation can be interpreted as a deviation from 
Gaussianity. Figure 7 shows the log-log plot of structure 
function in terms of time scaling (upper panel), expo- 
nents ^{q) (left lower panel) and C,{<l) (right lower panel) 
for the plasma fluctuations with / = 50mA. The present 
results are in agrement with our previous results derived 
that the plasma time series have multi- fractal nature [17]. 

The obtained expression for D(i)(A X.T) 
and I?(2)(Ax,r) (to avoid the overissue we just report 
the results of data for / = 50mA) are as follows 



TABLE 11. The values of moments, (a;"), and their errors 
for data set at different discharge current intensities. 





/„2\ 1 r,+5 


y -in+9 


\X ^ X lU 


50mA 


2.483 ± 0.001 


-0.395 ±0.417 


1.842 ± 0.062 


60mA 


2.415 ±0.001 


-3.330 ±0.123 


1.532 ±0.001 


100mA 


2.879 ± 0.001 


-17.500 ±1.990 


4.140 ±0.318 


120mA 


2.600 ± 0.001 


-3.140 ±0.179 


1.858 ± 0.008 


140mA 


2.095 ± 0.001 


-5.770 ±0.380 


1.431 ± 0.030 


180mA 


1.617 ±0.001 


-7.220 ± 2.850 


2.002 ± 0.765 


210mA 


1.383 ±0.001 


-0.458 ±0.139 


0.811 ±0.335 



D(i)(Aa;,r) = -(0.45 ± 0.03)Aa; 

D'^^^ (Ax, t) = (0.04 ± 0.01) Aa;2 (26) 

consequently, using Eqs. (26) and (24), the scaling expo- 
nent is determined as 

^{q) = (0.45 ± 0.03)g - (0.04 ± 0.01)q{q - 1) (27) 

As shown in the lower left panel of Figure 7, the above 
function for ^(g) (solid Hue) with the shaded area corre- 
sponds to 68.3% confidence interval derived by Marko- 
vian analysis of increments has an acceptable confidence 
level to experimental results (filled symbol). 

Finally we check the Gaussian nature of the PDFs of 
reconstructed and detrended time series. For a Gaussian 
distribution, all the even moments are related to the sec- 
ond moment through (cc^") = Jr^(a;^)" (e.g., for n = 2, 
(x^) = 3(x2)2), while the odd moments are zero iden- 
tically. We can directly check the relation between the 
higher moments for the plasma fluctuations data at dif- 
ferent value of discharge current intensities with second 
moment. The values of moments and their variances cal- 
culating directly from data are summarized in Table II. 

Let us examine the predictions for the moments of the 
plasma fluctuations via the Fokker-Planck equation, and 
compare their values with the direct evaluation repre- 
sented in Table II. Using the general Kramers-Moyal ex- 
pansion, Eq. (12), which is also valid for the probability 
density p{x,t), differential equations for the n-th order 
moments can be derived. By multiplication of the both 
side of Eq. (12) with x" and integration with respect to 
X, we can obtain evolution of different moments of data 
set as 

, OO (. + 00 / 3 \ '' 



= Er— Im/ x"-''D^''\x,t)p{x,t)dx 

k=l 

-E7;r^(-""'^*'H^,0) (28) 

k=l ^ '' 

We put n = 4 in the above equation and find the equation 
for the fourth moment as follows 

^{x^{t)) = A{D'^^\x)x^{t)) + 12{D^^\x)x^{t)) 
at 

+2A{D^^^x)x{t)) +2A{D'-^\x)). (29) 

The third and fourth Kramers-Moyal's coefficients for 
the data set are reported in Table III. We should point 
out that the values of \D^^^\ and \D^'^^\ are less than 
For the stationary case, all the moments of fluc- 
tuations are time independent and the left-hand side of 
Eq. (29) vanishes, so 

{x') = [a2(/)±a2(/)](x2)' + [a3(/)±a3(/)](x3)y^ 

(30) 

where a2{I) determines the coefficient of Kurtosis quan- 
tity and (/) shows its variance. Also a^il) determines 
the coefficient of Skewness and ^^{I) indicates its error. 
Skewness measures the asymmetry of probability den- 
sity function and Kurtosis determines the statistic of rare 
events in the processes. In generally they may depend to 
the discharge current intensity, /. Using the results rep- 
resented in Table III, for each case of fluctuations a„ (/) 
and its variance are given in Table IV. 

As we mentioned before, for exact Gaussian process, 
we should have 

"2(/) = ^=3.0 (31) 

«3(/) = 0.0 (32) 

If a2{I) > 3.0 means that probability density function 
has fat tail and rare events have more chance to occur. 
While for a2{I) < 3.0, the tails of probability density 
function is heavy than the Gaussian distribution. Ac- 
cording to the values of a2{I) and 03(1), we find that 
probability density function of data set is deviated from 
Gaussian. 

IV. SUMMARY AND CONCLUSION 

We have studied the stochastic nature of the electrical 
discharge current fluctuations in the Helium plasma as a 
working gas. Due to the sinusoidal trend we should dis- 
tinguish the intrinsic fluctuations from nonstationarity to 
infer valuable statistical properties of data set. Therefore 
for every stochastic analysis to avoid detection of spuri- 
ous statistical properties we had to first clean the data 



TABLE III. The values of third and fourth Kramers-Moyal 
coefficients for data set at different discharge current intensi- 
ties. 







L»W(a;) 


50mA 


_n 007 — oyn T 
-0.019 


n nnQ 4- ooq r -1- 01 

-0.001 a;^ 0.001 a;"* 


60mA 


—0 024- T — 002 
-0.010 


001 — 020 r -1- 01 
+0.001 + 0.003 x^ 


1 00m A 


—0 025 T — 001 r-^ 
-0.009 x^ 


002 + 001 r 4- 01 3 
+0.002 x^ 


120mA 


-0.013 X - 0.005 x^ 


0.001 +0.006 + 0.001 x-^ 


140mA 


-0.008 x - 0.003 x^ 


0.0009 + 0.003 x^ + 0.001 x'^ 


180mA 


-0.007 X - 0.003 


0.0007 + 0.002 x^ + 0.001 x^ 


210mA 


-0.009 X - 0.002 x^ 


0.0008 + 0.005 x^ 



set, then apply the statistical method to analysis the sig- 
nal. We have applied the Fourier-Detrended Fluctuations 
Analysis method to extract sinusoidal trend. 

We showed that how the mathematical framework of 
Markov processes can be applied to develop a success- 
ful statistical description of the plasma fluctuations. We 
have analyzed detrended data via Markovian method. 
The Markov time scale, as the characteristic time scale of 
the Markov properties of the electrical discharge current 
fluctuations, was obtained. According to the theory of 
the stochastic process, the electrical discharge current at 
time scales larger than the Markov time scale can be con- 
sidered as a Markov process. This means that the data 
located at the separations larger than the Markov time 
scale can be described as a Markov chain. It is found 
that Markov time scale, tMarkov increases by increasing 
the current intensity in the plasma. This means that 
the memory of charged particles in the plasma increase 
as current intensity increases. It is due to the fact that 
particles become more energetic, therefore they can pene- 
trate deeper in the plasma without considerable deviation 
from the initial way. Then every new situation keeps the 
memory of some nearby conditioning corresponds to its 



TABLE IV. The values of a's coefRcients and their vari- 
ances, cr's, for data set at different discharge current intensi- 
ties. 





02 


0"2 


as 


0-3 


50mA 


3.14 


0.38 


0.18 


0.21 


60mA 


3.81 


0.70 


-1.17 


0.51 


lOOrnA 


2.86 


0.53 


-3.68 


0.66 


120mA 


2.77 


0.98 


2.25 


0.96 


140mA 


3.05 


0.56 


-2.09 


0.35 


180mA 


3.67 


0.85 


-0.01 


0.28 


210mA 


3.08 


0.95 


-4.18 


0.95 



Markov time scale. In other words by increasing the cur- 
rent intensity, effective scattering cross-section decreases 
so mean free time increases [53] as weh as corresponding 
Markov time scale. According to the Markovian nature 
of fluctuations, we demonstrated that, the probability 
density function of fluctuations satisfies a Fokker-Planck 
equation. The Langcvin equation, governing the evolu- 
tion of current fluctuations also has been given. 

To check the multifractality nature of the time series, 
we used the scaling properties of structure function and 
concept of Extended-Self-Similarity [54,55]. By comput- 
ing the corresponding Kramers-Moyal's coefficients, we 
find a good consistency between multifractality nature of 
original time series and the obtained exponents from the 
reconstructed Fokker-Planck equation. 

Acknowledgements Authors would like to thank S. 
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